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Abstract 

The nuclear physics dependence of the exotic (u~,e~) conversion branching ratio i? Me - 
for the experimentally most interesting nuclei 208 Pb and 48 Ti, is investigated in various 
nuclear models. The results thus obtained are combined with the new experimental limits 
extracted at PSI for these nuclei to put bounds on the elementary particle parameters 
entering R^ e - such as intermediate neutrino masses and mixing angles as well as relevant 
parameters of intermediate super symmetric particles (masses and mixing of s-fermions 
and neutralinos). 
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There is a plethora of processes predicted by many extensions of the standard model like 
grand unification, supersymmetric theories etc., which violate the lepton family quantum num- 
bers [0-0. Among them the neutrinoless conversion of a bound muon to an electron, 

(A,Z) + pi - e" + (A,Z)* (1) 

which violates the muon and electron numbers, stands out as one of the most prominent pro- 
cesses to search for lepton flavour non-conservation if it exists [f|-[f|. The recent experiments 
at PSI |5], |5| and TRIUMF [|7|], have up to now put only upper limit on the branching ratio 
R pe - = r(/i~, e _ )/r(/i~, i/p) in the values 

Rll- < 4.3 x KT 12 , for 48 Ti target § 

R™- < 4.6 x KT 11 , for 208 Pb target § 

We mention that, the limit set this year at PSI by using 208 Pb as target ||, improved over 
the previous limit |7j by an order of magnitude. By progressively improving the experimental 
sensitivity, it is expected in the next few years these limits to be pushed down by two to three 
orders of magnitude 0-0. 

On the theoretical side, up to now the elementary particle physics aspects of the (/i~,e~) 
conversion amplitudes have been investigated in the context of several models allowing lepton 
flavour violation [l|]-[f|, 0, O n the other hand, the nuclear form factors, which describe the 
rates of process (1) have been computed employing various nuclear models [n|-[IE|. In recent 



years some refinements of the quasi-particle RPA approximation, which have been popular in 
evaluating successfully the nuclear transition matrix elements entering processes like double 
beta decay (see e.g. ref. [jnj), double charge exchange |T3 etc, have also been employed in the 



investigation of (/i~,e~) conversion in the presence of nuclei [12 . 

In the present work, the p-p and n-n quasi-particle RPA is used in evaluating the nuclear 
matrix elements entering R pe - for the targets 208 Pb and 48 Ti, for various mechanisms leading to 



reaction (1) utilizing the method developed in ref. ||12j| . Our goal is, by combining the nuclear 
calculations with the experimental limits quoted above, to set constraints on the elementary 
particle parameters entering various models involving intermediate massive neutrinos or other 
exotic supersymmetric particles (s-leptons and neutralinos) as it will be discussed below. We 
perform explicit calculations of the nuclear matrix elements of Pb and Ti in the appropriate 
values of the momentum transfer by taking into account the muon binding energy entering 
both the phase space and the nuclear transition form factors. Furthermore, we calculate the 
contributions coming from Z-exchange diagrams which have not been included in previous 
studies [O, O. We have paid special attention to the nucleus 208 Pb, not only because it is 



currently used as target at PSI but, as we have previously shown |T3[, the rate of the (//~ , e~) 
conversion appears to attain a maximum in this region. We mention that, since 208 Pb is a 
doubly closed shell nucleus, it needs a special treatment in the context of quasi-particle RPA 

The nuclear calculations performed in this work rely on the effective interaction (/x , e ) 
conversion Hamiltonian constructed in the framework of common extensions of the standard 
model involving mixing of intermediate neutrinos - gauge bosons or neutralinos - s-leptons 
(supersymmetry). Below we describe in brief the essential formalism required for our study. 

1. Intermediate neutrinos - gauge bosons. 

It is well known that, within the minimal standard model, the strengths of the flavour 
changing interactions are related to the neutrino masses. In this spirit, one writes the weak 
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neutrino eigenstates u e , u^, ... in terms of the mass eigenstates Uj (light neutrinos) and Nj (heavy 
neutrinos) as 

n n 

E^ + E^- ( 2 ) 



V e = 

3=1 3=1 



E^' + E^i (3) 

J=l J=l 

where n=number of generations and l/j^j are the elements of the charged lepton current 
mixing matrix for light neutrinos Uj associated with the electron and the muon (see figs. 1,2). 
The corresponding quantities with superscript (2) refer to heavy neutrinos Nj. 

In order to deduce the effective interaction Hamiltonian of the (n~ , e~) process at the nuclear 
level needed for our study, one starts from the weak vector and axial vector quark currents, 

j a = E^Ta(1 -j 5 )qj, Jx = E^TaCI-Ts)^- (4) 



(J\ neutral, J x charged currents) where qj,q'j = u,d quarks (we neglect strange quark contri- 
bution into the nucleon and consider only left-handed currents). 

The next step is, to write at the nucleon level the hadronic current for each specific mecha- 
nism by taking the matrix elements of eq. (f|) using the appropriate quark model wave function 
for the nucleon. In this work we assumed a non-relativistic nucleon wave function and studied 
the following types of mechanisms: 

i) Photonic diagrams: The (fi~,e~) conversion can proceed via the diagrams of figs. 1(a) 
and 2(a), i.e. by exchange of a virtual photon, provided the upper vertex is allowed. The 
hadronic vertex is the usual electromagnetic coupling. At the nucleon level it can be written as 

J? = N plx N p = iV 7A i (1 + r 3 ) N (photonic) (5) 

(N p , N represent the proton, nucleon spinors) This current involves equal isoscalar and isovector 
components. 

ii) Non-photonic diagrams: There is a plethora of such diagrams. The most obvious are the 
one which is mediated by Z-particle exchange shown in figs. l(a),(b) and 2(a),(b), and the box 
diagrams involving intermediate massive neutrinos and W-bosons as in fig. 1(c) or mediated 
by s-leptons and neutralinos as in fig. 2(c). In such cases the hadronic current takes the form 

J? = N lx \ [(3 + f v (3r 3 ) - (f v f3" + f A P'r 3 ) l5 \ N (nan - photonic) (6) 

The parameter takes the value unity except in Z-exchange when it is zero. For the box 
diagrams we have (3 = (3 . In the case of W-boson exchange we have (3 = 5/6 while for 
intermediate s-leptons we get (3 = 0.6. For Z-exchange we find 

f 3 ' = v^w = 6 - 90 W 

13 = 3(1 ^— ) = -3.46 (8) 

H v 2sin 2 6 w ' v ; 

In the above equations (3 = Pi/ Po with p (Pi) being the isoscalar (isovector) quantities at the 
quark level, /y, /U are the vector, axial vector static nucleon form factors (Ja/ fv = 1.24). 
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The effective Lagrangian at the nucleon level, which involves both photonic and non- 
photonic contributions can be cast in the form 

M = ^,4" + ^(Vf 0) 

where a is the fine structure coupling constant, ( is given by 

G F m 2 

( = — -j=J- (W — boson exchange) (10) 

and q is the momentum transfer which in a good approximation is written as 

q = mfl -e-(E f -E gs ) (11) 

with e the muon binding energy and Ef (E gs ) the energy of the final (ground) state of the 
nucleus. 

In eq. and jh^ represent the leptonic currents corresponding to each mechanism i.e. 

i) photonic mechanism: 

jfa = u(Pe) [ {f Ml + Tb/iSl)^— + (fEO + 75/Mo) lu (d* ~ ^f) ]u(p M ) (12) 

ii) Non-photonic mechanism: 

J(2) = u{Pe)l X \ (fl + MsMPm) (13) 

where p e , p^ the lepton momenta and /eq, /ei, fhio, /mi, /i, /2 parameters depending on the 
assumed gauge model. 

In the last step, the effective interaction Hamiltonian, which converts a muon to an electron 
in the nucleus (A, Z) , is obtained by summing over all type of eq. (§) single nucleon contri- 
butions. Behind this summation lies the assumption that the A nucleons interact individually 
with the muon field (impulse approximation). 

2. Intermediate neutralinos - s-leptons (supersymmetric model). 

The supersymmetry (SUSY) associates to each particle its superpartner (s-particle). In this 
way, the fundamental particles of nature are essentially doubled and as a consequence several 
new mechanisms, some of which violate the lepton flavour conservation, are predicted. Some 
typical supersymmetric diagrams leading to the (fi~,e~) conversion are shown in fig. 2. 

In the effective Hamiltonian of eq. the hadronic current is given by eq. @ by putting 
the values of (3 discussed before. For the leptonic currents, in the simplest version of this model, 
the photonic and non-photonic amplitudes are simply related |10[|. One finds: 



(47ra)/ M i = -(47ra)/ E1 = -— / (14) 

(4na)f E0 = -(47ra)/ M o = ~f (15) 

fi = -h = ^Pof (16) 

/ = (17) 
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fj= { *^ (18) 
m 2 

C = 1.0 (19) 

where m 2 is the average of the square of the masses of the s-fermions entering the loop and 
/?o = 5/9. The parameter fj depends on the details of the model [ITU . 



The Z-exchange contribution in SUSY models like the one just discussed yields [pOfl 
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fi = ~h = « 2 ^^ (20) 

where the parameter £ depends on the details of the model. It vanishes if one ignores the 
Higgsino components of the neutralinos as e.g. in the case of pure photino. It also vanishes if 
the probabilities of finding the two Higgsinos in the neutralino are equal. In the model #2 of 
ref. |TJ| it takes the value £ = 2.8 x 10~ 2 . Clearly, the Z-exchange is more favored in models in 
which 7h 2 is larger than m%. Accurate formulas and details will be presented elsewhere |2Q| . 

Experimentally, the branching ratio R^- is quite interesting quantity. From a theoretical 
point of view, in general it is not easy to formulate an expression containing contributions 
from both mechanisms, photonic and non-photonic ones. For this reason, the photonic and 
non-photonic contributions are usually discussed separately. 

In the case of the coherent channel, however, one can separate the dependence on the nuclear 
physics from the leptonic form factors in the branching ratio i? Me -, provided that the nuclear 
vertex is calculated in the Born approximation and the density of the muon bound state inside 
the nucleus is approximated by a mean wave function (see below). Under these assumptions 
the muon capture rates can well be described by the Goulard-Primakoff function fcp(A,Z) 
[ PH . In general, it is impossible to separate the nuclear structure aspects from the elementary 
particle parameters. This can, however, approximately be done in the models discussed above 
[[Toll . Thus one can write 

R»e-=P1 (21) 

The function j(A, Z) contains all the nuclear dependence of R^ e - and the quantity p is inde- 
pendent on nuclear physics. They are defined as follows: 

i) Neutrino mixing models: 

E ePe \ME\ 2 . . 

_ _eFe 1 1 ^ 

1 ml G 2 Zf GP (A,Z) 1 1 

(G 2 ~ 6, p e = — q) where \ME\ represent the nuclear matrix elements given in terms of 
the proton (Fz(q 2 )) and neutron (F/v(g 2 )) nuclear form factors. In the special case of photonic 
diagrams \ME\ is proportional to the elastic or inelastic form factor Fz(q 2 ) depending on which 
final state is populated. Then, 7 takes the form 

EePe Z\F z (q 2 )\ 2 
lph ml G 2 f GP (A,Z) [Z6) 

The quantity p in eq. fl2Tj), in the case of left-handed currents only, takes the form 



1 a \2 I /mi + /sol 2 + |/bi + /mo I 
p = (47ra) 



2^2 



— — I — i^Vu + Vn\ 2 (photonic diagrams) (24) 
647i 2 77% 
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|/3o/i| 2 +|/3o/ 2 



m 



2 



= — 4 (Gj?m f ) 1 30 — ^-f] u + r] N \ (non — photonic diagrams) (25) 

The lepton violating parameters associated with intermediate light (j] u ) or heavy (t]n) neutrinos 
depend on the gauge model considered. They are given by 

*=?w^. ^=E^ , ^- , *5(" 2i ' i 3 +3 ) (26) 



ii) Neutralino and s-lepton mixing models: Flavour violation in these models occurs via the s- 
lepton mixing provided that it is different from the corresponding charged lepton mixing. It also 
involves intermediate neutralinos. Following jnj we will assume that the dominant contribution 
comes from the photino (SUSY partner of the photon). Under some plausible assumptions Q 
and neglecting the Z-exchange we find that the function j(A, Z) can be written as 

,13 1NF N ^ 
The quantity p in this model can be written as 

a 4 



P 288 (G F fh 2 ) 2 ^ ^ 



where rj is given by eq. fljjD - In the case of Z-exchange |20| we find 



7 = 0.053(1 - 14.04 y^ffl* (29) 



The quantity p is now given by 



z 



p = 5.5 x 10~ 4 — — - \r]\ 2 (30) 



(Gftu 



We stress that for the coherent process, in the models discussed above, the only variable of 
the elementary particle sector, which enters the function 7 is the parameter j3 = f3\/ (3q. Once 



7(A, Z) is known, e.g. by nuclear model calculations, from eq. (^T|) one can extract information 
about the interesting parameter p from the experimental data and compare it with the value 
predicted for various mechanisms by the gauge models. Such kind of calculations we include 
below. 

The task of studying the nuclear physics aspects of the exotic (/i~,e~) conversion process 
is to evaluate the function j(A,Z) of eqs. (p2|), (p3[), (|27|) and (p9|). In general, the nuclear 



matrix elements entering "f(A, Z) depend on the final nuclear state populated during process 
(1). What, however, is more interesting is the coherent process, since then the rate is free of 
the background from bound muon decay p. Furthermore, extensive calculations have shown 



that the coherent channel dominates the (pT ', e~) conversion throughout the periodic table [13 
(for the isotopes studied in the present work see table 2 below). 

If we assume an average value for the muon wave function, so that it cancels in the branching 
ratio, the coherent matrix elements needed in the computation of j(A, Z) can be written as 
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M 2 coh (q 2 



Z 2 F^) (31) 



where Q(/3), for the models discussed in the present work, takes the values 



W) 



3 + fvP 



0, photonic diagrams, (3 = 3 

13/23, W — boson exchange, (3 = 5/6 , . 

2/3, SUSY Photonic + Box diagrams, (3 = 3/5 

-14.04, SUSY Z - exchange, (3 = -3.46 



Notice that, in the photonic case M 2 oh = Z 2 F^(q 2 ) (only protons of the considered nucleus 
contribute, see eq. fl5|)). 

In table 1, we show the results of the nuclear form factors and the coherent matrix elements 
for 208 Pb and 48 Ti, calculated in the context of quasi-particle RPA using the method of ref. 
[O]. For comparison, in this table the results of ref. [[□]] (shell model) and those of ref. [13 



(local density approximation, LDA) are also presented. The quasi-particle RPA form factors 



are in good agreement with experimental electron scattering data ]2B| also listed in table 1. 

For the benefit of the reader we provide here a brief description of the main steps followed 
in the computational procedure. At first, we have chosen an appropriate model space: for 208 Pb 
18 levels above the core 59° Sn and for 48 Ti 16 levels without core (for the calculations of 48 Tz 
in ref. |T2] we used a smaller model space, consisted of only 10 levels without core). In order to 



satisfy the convergence of BCS equations in the case of the doubly closed shell nucleus 208 Pb, we 
determined the strength parameters of pairing for protons and neutrons from the experimental 
proton and neutron pairing gaps of the (N — 2, Z + 2) neighbor nucleus 124-P0, following the 
procedure used recently in the study of double beta decay |L7|. 

We have also evaluated the incoherent rate, i.e. the matrix elements from the initial (ground) 
state (a + state) to every excited state \f > included in the chosen model space of the studied 
nucleus, by calculating explicitly the contribution of each individual channel, for photonic 
and non-photonic diagrams. These results are denoted as Mf nc and are shown in table 2 
separately for the vector (My) and axial vector (Ma) components of the hadronic current eq. 
(^|) (M 2 nc = My + Ma). Obviously, in the photonic mechanism only the vector component 
contributes, but in the non-photonic one both the vector and axial vector components give 
non-zero contributions. For each isotope studied we found that, the main contribution to 
the incoherent rate comes from the low-lying excited states and that high-lying excited states 
contribute negligible amounts. Such state-by-state calculations of the incoherent channel for a 



set of nuclei throughout the periodic table are discussed in ref. |T2|, ^2 . 

From the coherent and total matrix elements, M 2 ot = M 2 oh + Mf nc , another useful quantity 
of the (/i~,e~) conversion, the ratio n = M 2 oh /M 2 ot which expresses the portion exhausted by 
the coherent rate in the total branching ratio R^e-, can be obtained (see table 2). We found 
that, for both isotopes and all mechanisms studied, n > 90%, which implies that the coherent 
rate dominates the (/i~, e~) conversion process. In earlier calculations the ratio n was estimated 
FJ|] to be i] « 83 % in Cu region and smaller in lead region while in ref. [13] we found rj > 90% 



throughout the periodic table. Our present results agree well with those of ref. |13j . 

We should note that, in the present study the muon binding energy e& in eq. ([LI]) has been 
taken into consideration. This property affects significantly the nuclear form factors, especially 
for heavy nuclei like 208 Pb. In light nuclei e& is negligible (for Ti e& = 1.450), but in heavy 
elements it becomes significant (for Pb e& = 10.450). In addition, the factor E f /p e jm l L in eq. 
which takes into account the phase space in the transition matrix elements, depends on 
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the muon binding energy. By ignoring e^, this factor becomes unity. For 48 Ti this factor is 
equal to 0.97 and the neglect of e& is a good approximation. For 208 Pb, however, this factor is 
equal to 0.81 which means that, in heavy nuclei the dependence on e& of the phase space cannot 
be ignored. 

The values for 7 (A, Z) calculated from the coherent QRPA matrix elements are shown in 
table 3 and compared with the results of refs. [11, 13| . We mention that, the variation of 



7(^4, Z) through the periodic table studied in ref. [JTTJ] exhibits a strong dependence on the 
neutron excess (N — Z) which mainly reflects the dependence on (N, Z) of the total muon 
capture rate. 

By putting the experimental limits |5|, || in eq. (|21~D we determine upper bounds on the 
parameter p for photonic, Z-exchange and W-box diagrams both in conventional extensions of 
the standard model as well as SUSY theories (see table 3). We note that the values of p for 
i8 Ti are slightly improved over those of table 2 of ref. ]l!|, but those of 208 Pb are appreciably 



smaller than those of sect. 5.1 of ref. [11]. This big difference is due to the fact that, in the 
present work we have used the experimental limit of ref. || . This experiment at PSI improved 
on the previous limit of 208 Pb 0, which was used in ref. [ TT[] , by an order of magnitude. For 



both nuclei we found that, the limits obtained from various nuclear models do not significantly 
differ from each other (we should stress that, the different limits for 208 Pb in the shell model 



results of table 3, are due to the neglect of e& in eq. (11) when calculating the nuclear form 
factors as mentioned above; its consideration gives similar results to those of LDA and QRPA 
for both mechanisms). This implies that all nuclear models studied here give about same values 
for p. 

One can use the limits of p to parametrize the muon number violating quantities entering 
eqs. (|24l) , (|25D , (|28"D and (|30"D and also to compare them directly to the value given from gauge 
models. As an example, we quote the value of p = 8.2 x 10~ 18 obtained in the supersymmetric 
model of ref. |TTJ discussed above. This prediction of p is considerably smaller compared to 
the values listed in table 3, which were extracted from experiment. 

In summary, in the present letter we have studied the nuclear physics part of the branching 
ratio of the exotic p~ — > e~ conversion R^ e - in 208 Pb and 48 Ti nuclei. These two isotopes are 
the most interesting nuclear targets to search for lepton flavour violation. 208 Pb is currently 
used at PSI in the SINDRUM II experiment. 

We have calculated in a reliable way the appropriate nuclear matrix elements entering the 
branching ratio R^ e - in the context of the quasi-particle RPA and compared them with the 
results given from other nuclear models. We found that, the coherent (p~,e~) rate, which is 
measured from experiments, dominates the branching ratio R^ e - for both isotopes but it is 
more pronounced in the heavy nucleus 208 Pb. 

From the calculations of the nuclear part of the branching ratio R^ e - , using the new ex- 
perimental limits, especially those for 208 Pb, we were able to estimate upper limits for the 
elementary sector part of the pT — > e~ rate in common extensions of the standard model, 
which allow lepton flavour violation. These limits are very useful to fix the lepton violating 
parameters and test the various gauge models. 
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Table 1. Quasi-particle RPA results for proton, neutron nuclear form factors and coherent 
nuclear matrix elements squared for the four cases discussed in the text: photonic and non 
photonic W-boson exchange as well as the two SUSY cases (photonic + box and Z-exchange). 
For comparison the experimental charge form factors {F^ xp ) P3| and the calculations of shell 



model [ED] and LDA [[13] are shown. Note that in the Z-exchange the factor 0.053 of eq. (p9[) 



has not been included. 





Shell Model 


LDA 


Exper. 




Quasi-particle RPA Results 




(A,Z) 


F z 


F N 


Fz 


F N 


pexp 


F z 


F N M p 2 ft M^_ ex M 2 _ lep 


M z - ex 


i8 Ti 


.543 


.528 


.528 


.506 


.532 


.537 


.514 139.6 375.2 429.6 


30918.0 


208 pfo 


.194 


.139 


.250 


.220 


.242 


.271 


.214 494.7 1405.2 1618.8 


127214.1 



Table 2. The (/i~, e~) conversion matrix elements (coherent, incoherent, total) and the ratio 
7] = M 2 oh /M 2 ot for photonic and W-exchange (non-photonic) mechanisms involving intermediate 
neutrino mixing. Results for the experimentally most interesting nuclei, 208 Pb and 48 Ti, are 
presented. 





Photonic fi — > e 


Mechanism 


W-exchange \i — > 


e Mechanism 




M 2 coh 


M 2 

inc 


M tot V(%) 


M 2 coh 


M 2 


M tot V(%) 


(A,Z) 




My M A 






My M A 




48 T , 


139.6 


14.1 - 


153.7 90.8 


375.2 


10.7 2.6 


388.5 96.6 


208p^ 


494.7 


14.0 - 


508.7 97.2 


1405.2 


19.0 6.8 


1431.0 98.2 
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Table 3. The new limits on the elementary sector part of the exotic /x — e conversion 
branching ratio extracted by using eq. Q2TJ) and the recent experimental data for the nuclear 
targets 208 Pb, iS Ti 
7 (A, Z) of eqs. 



H P • The nuclear part of the branching ratio, described by the function 
Q23D, Q27D and is also shown. 









48 T , 






208p^ 




Method 


Mechanism 


l(A, Z) 




r 




l(A, Z) 


o 

r 






Photonic 


9.42 


< 4.6 xlO" 


-13 


17.33 


< 2.7 xlO" 


-12 


QRPA 


W-boson exchange 


25.31 


< 1.7 xlO" 


-13 


49.22 


< 0.9 xlO' 


-12 




SUSY s-leptons 


25.61 


< 1.7 xlO" 


-13 


49.50 


< 0.9 xlO' 


-12 




SUSY Z-exchange 


110.57 


< 0.4 xlO" 


-13 


236.19 


< 0.2 xlO" 


-12 


LDA 


Photonic 


9.99 


< 4.3 xlO" 


-13 


17.84 


< 2.6 xlO" 


-12 




W-boson exchange 


26.60 


< 1.6 xlO" 


-13 


55.53 


< 0.8 xlO" 


-12 


SM 


Photonic 


9.74 


< 4.4 xlO" 


-13 


10.42 
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Figure Captions 

Figure 1. Typical diagrams entering the neutrinoless (/J~,e~) conversion: photonic 1(a), 
Z-exchange l(a),(b) and W-boson exchange 1(c), for the specific mechanism involving interme- 
diate neutrinos. 



Figure 2. SUSY diagrams leading to the (/i~,e - ) conversion: photonic 2(a), Z-exchange 
2(a),(b) and box diagrams 2(c), in a supersymmetric model with charged s-lepton and neutralino 
mixing are shown. Note that, the Z-exchange in 2(b) as well as in fig. 1(b), comes out of 
electrically neutral particles (photon exchange does not occur in these diagrams). These Z- 
exchange diagrams may be important (those of figs. 1(a) and 2(a) are suppressed by m^/ml). 
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